Charge correlations and optical conductivity in weakly doped antiferromagnets

We investigate the dynamical charge-charge correlation function and the optical conductivity in weakly doped antiferromagnets using Mori-Zwanzig projection technique. The system is described by the two-dimensional t-J model. The arising matrix elements are evaluated within a cumulant formalism which was recently applied to investigate magnetic properties of weakly doped antiferromagnets. Within the present approach the ground state consists of non-interacting hole quasiparticles. Our spectra agree well with numerical results calculated via exact diagonalization techniques. The method we employ enables us to explain the features present in the correlation functions. We conclude that the charge dynamics at weak doping is governed by transitions between excited states of spin-bag quasiparticles.Comment: 5 pages, 2 figures, to appear in Europhys. Letter

Pairing of 1-hexyl-3-methylimidazolium and tetrafluoroborate ions in n-pentanol

Molecular dynamics simulations are obtained and analyzed to study pairing of 1-hexyl-3-methylimidazolium and tetrafluoroborate ions in n-pentanol, in particular by evaluating the potential-of-mean-force between counter ions. The present molecular model and simulation accurately predicts the dissociation constant Kd in comparison to experiment, and thus the behavior and magnitudes for the ion-pair pmf at molecular distances, even though the dielectric constant of the simulated solvent differs from the experimental value by about 30%. A naive dielectric model does not capture molecule structural effects such as multiple conformations and binding geometries of the Hmim+ and BF4- ion-pairs. Mobilities identify multiple time-scale effects in the autocorrelation of the random forces on the ions, and specifically a slow, exponential time-decay of those long-ranged forces associated here with dielectric friction effects.Comment: 5 pages, 7 figures. V2: Figs. 4 & 7 redrawn for better visual clarity with log-scales. No change in results. In press J. Chem. Phys. 201

Exact dynamics of dissipative electronic systems and quantum transport: Hierarchical equations of motion approach

A quantum dissipation theory is formulated in terms of hierarchically coupled equations of motion for an arbitrary electronic system coupled with grand canonical Fermion bath ensembles. The theoretical construction starts with the second--quantization influence functional in path integral formalism, in which the Fermion creation and annihilation operators are represented by Grassmann variables. Time--derivatives on influence functionals are then performed in a hierarchical manner, on the basis of calculus--on--path--integral algorithm. Both the multiple--frequency--dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting formalism is in principle exact, applicable to interacting systems, with arbitrary time-dependent external fields. It renders an exact tool to evaluate various transient and stationary quantum transport properties of many-electron systems. At the second--tier truncation level the present theory recovers the real--time diagrammatic formalism developed by Sch\"{o}n and coworkers. For a single-particle system, the hierarchical formalism terminates at the second tier exactly, and the Landuer--B\"{u}ttiker's transport current expression is readily recovered.Comment: The new versio

Forcing anomalous scaling on demographic fluctuations

We discuss the conditions under which a population of anomalously diffusing individuals can be characterized by demographic fluctuations that are anomalously scaling themselves. Two examples are provided in the case of individuals migrating by Gaussian diffusion, and by a sequence of L\'evy flights.Comment: 5 pages 2 figure

Collective motion of binary self-propelled particle mixtures

In this study, we investigate the phenomenon of collective motion in binary mixtures of self-propelled particles. We consider two particle species, each of which consisting of pointlike objects that propel with a velocity of constant magnitude. Within each species, the particles try to achieve polar alignment of their velocity vectors, whereas we analyze the cases of preferred polar, antiparallel, as well as perpendicular alignment between particles of different species. Our focus is on the effect that the interplay between the two species has on the threshold densities for the onset of collective motion and on the nature of the solutions above onset. For this purpose, we start from suitable Langevin equations in the particle picture, from which we derive mean field equations of the Fokker-Planck type and finally macroscopic continuum field equations. We perform particle simulations of the Langevin equations, linear stability analyses of the Fokker-Planck and macroscopic continuum equations, and we numerically solve the Fokker-Planck equations. Both, spatially homogeneous and inhomogeneous solutions are investigated, where the latter correspond to stripe-like flocks of collectively moving particles. In general, the interaction between the two species reduces the threshold density for the onset of collective motion of each species. However, this interaction also reduces the spatial organization in the stripe-like flocks. The most interesting behavior is found for the case of preferred perpendicular alignment between different species. There, a competition between polar and truly nematic orientational ordering of the velocity vectors takes place within each particle species. Finally, depending on the alignment rule for particles of different species and within certain ranges of particle densities, identical and inverted spatial density profiles can be found for the two particle species.Comment: 16 pages, 10 figure

Nuclear quantum effects in solids using a colored-noise thermostat

We present a method, based on a non-Markovian Langevin equation, to include quantum corrections to the classical dynamics of ions in a quasi-harmonic system. By properly fitting the correlation function of the noise, one can vary the fluctuations in positions and momenta as a function of the vibrational frequency, and fit them so as to reproduce the quantum-mechanical behavior, with minimal a priori knowledge of the details of the system. We discuss the application of the thermostat to diamond and to ice Ih. We find that results in agreement with path-integral molecular dynamics can be obtained using only a fraction of the computational effort.Comment: submitted for publicatio

Lower bounds for the conductivities of correlated quantum systems

We show how one can obtain a lower bound for the electrical, spin or heat conductivity of correlated quantum systems described by Hamiltonians of the form H = H0 + g H1. Here H0 is an interacting Hamiltonian characterized by conservation laws which lead to an infinite conductivity for g=0. The small perturbation g H1, however, renders the conductivity finite at finite temperatures. For example, H0 could be a continuum field theory, where momentum is conserved, or an integrable one-dimensional model while H1 might describe the effects of weak disorder. In the limit g to 0, we derive lower bounds for the relevant conductivities and show how they can be improved systematically using the memory matrix formalism. Furthermore, we discuss various applications and investigate under what conditions our lower bound may become exact.Comment: Title changed; 9 pages, 2 figure

Accelerating the convergence of path integral dynamics with a generalized Langevin equation

The quantum nature of nuclei plays an important role in the accurate modelling of light atoms such as hydrogen, but it is often neglected in simulations due to the high computational overhead involved. It has recently been shown that zero-point energy effects can be included comparatively cheaply in simulations of harmonic and quasi-harmonic systems by augmenting classical molecular dynamics with a generalized Langevin equation (GLE). Here we describe how a similar approach can be used to accelerate the convergence of path integral (PI) molecular dynamics to the exact quantum mechanical result in more strongly anharmonic systems exhibiting both zero point energy and tunnelling effects. The resulting PI-GLE method is illustrated with applications to a double-well tunnelling problem and to liquid water

Specific heat anomalies of open quantum systems

The evaluation of the specific heat of an open, damped quantum system is a subtle issue. One possible route is based on the thermodynamic partition function which is the ratio of the partition functions of system plus bath and of the bath alone. For the free damped particle it has been shown, however, that the ensuing specific heat may become negative for appropriately chosen environments. Being an open system this quantity then naturally must be interpreted as the change of the specific heat obtained as the difference between the specific heat of the heat bath coupled to the system degrees of freedom and the specific heat of the bath alone. While this difference may become negative, the involved specific heats themselves are always positive; thus, the known thermodynamic stability criteria are perfectly guaranteed. For a damped quantum harmonic oscillator, instead of negative values, under appropriate conditions one can observe a dip in the difference of specific heats as a function of temperature. Stylized minimal models containing a single oscillator heat bath are employed to elucidate the occurrence of the anomalous temperature dependence of the corresponding specific heat values. Moreover, we comment on the consequences for the interpretation of the density of states based on the thermal partitionfunction.Comment: 7 pages, 6 figures, new title and some modifications of the main tex

Efficient stochastic thermostatting of path integral molecular dynamics

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high-frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently-developed colored-noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nos\'e-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.Comment: Accepted for publication on JC